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Abstract - We study a model of Tensorial Group Field Theory (TGFT) on R 3 from the point of 
view of the Functional Renormalisation Group. This is the first attempt to apply a renormalisation 
procedure to a TGFT model defined over a non-compact group manifold. IR divergences (with 
respect to the metric on R) coming from the non-compactness of the group are regularised via 
compactification, and a thermodynamic limit is then taken. We identify then IR and UV fixed 
points of the RG flow and find strong hints of a phase transition of the TGFT system from a 
symmetric to a broken or condensate phase in the IR. 


Introduction. — Group Field Theories [1 {5J (GFTs) 
are a particular class of quantum field theories with fields 
defined over a group manifold and characterised by combi- 
natorially non-local interaction terms. This combinatorial 
non-locality makes the Feynman diagrams of the theory 
stranded diagrams dual to cellular complexes (simplicial 
complexes in the simplest constructions) [l,3j. GFT’s his¬ 
torically were born from an attempt to generalise matrix 
models [6] of 2d gravity to higher dimensions in the form 
of tensor models 7■ 9|. These models were soon enriched 
with group theoretic data in such a way that the Feynman 


sentation of the same amplitudes [16] . 

One key open issue of all of the above quantum grav¬ 
ity approaches, and GFTs in particular, is the emergence 
of a continuous geometry out of the discrete and quan¬ 
tum pre-geometric structures defining the formalisms, and 
of General relativity as an effective description of their 
(collective) dynamics in the same approximation. The 
study of the GFT collective dynamics and of the asso¬ 
ciated continuum limit is therefore crucial. Moreover, one 
suggested scenario for the emergence of spacetime and ge¬ 
ometry out of such quantum gravity models involves a 



models) endowed with a discrete geometric interpretation 
(which requires appropriate group theoretic data and suit¬ 
able choices of dynamics) there is also a direct link with 
simplicial quantum gravity path integrals, first manifested 
in the semiclassical analysis of GFT Feynman (spin foam) 
amplitudes, where one recovers the Regge action [15], and 
then shown to be generically manifest in the flux repre¬ 


that it offers the possibility to address this issue with tools 
coming from standard quantum field theory. 

And indeed, an important line of recent developments 
has concerned the renormalisation of GFT models, since 
the renormalisation group is indeed the key tool to address 
both the quantum consistency of field theory dynamics 
and the definition of the continuum limit, aimed at a pre- 
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cise mapping of the phase diagram of the theory. Further¬ 
more, GFT renormalization is also one of the two main 
strategies to define and study the renormalization of spin 
foam models, the other being through a generalised lattice 
gauge theory approach [22]. Most work in this direction 
has concerned a particular class of GFTs, called Tensorial 
Group Field Theories (TGFT’s) (23 -33 , which incorpo¬ 


rate recent advances in the statistical analysis of colored 
tensor models 34 - 37 . In particular, in TGFT framework, 
fields are endowed with tensorial transformation proper¬ 
ties under the action of the group itself. The perturbative 
analysis of these field theories has been undertaken and 
a large set of models prove to be perturbatively renor- 
malisable and asymptotically free (see references above). 
However, understanding the continuum limit, including 
the phase diagram and phase transitions of the same mod¬ 
els requires the study of their non-perturbative properties. 

This being the goal, the Functional Renormalisation 


Group is an efficient framework to reach it 38 -41 


The FRG approach has been applied first to matrix 
models 42 - 44] (with the double scaling critical point re¬ 
interpreted as a fixed point in the RG flow). More recently, 
the FRG framework has been adapted and applied for the 
first time to TGFTs in |45l. The authors of |45l studied 


a rank-3 TGFT defined over a compact 17(1) group mani¬ 
fold. The /3-functions define a non-autonomous system in 
the cut-off N. Then, the authors studied two regimes of 
the cut-off, the large and small N (and also an intermedi¬ 
ate regime at fixed TV), where a proper notion of dimension 
of the couplings can be defined and an autonomous sys¬ 
tem of RG equations is obtained. The notion of UV or 
IR “fixed points” is then only loosely (i.e. asymptotically) 
defined, as the existence of a trajectory from a UV to a 
IR fixed point becomes more difficult to ascertain. This is 
not surprising nor problematic per se, and it simply sig¬ 
nals the presence of an additional scale in the formalism, 
here the size of the group manifold on which the fields are 
defined. In fact, the same feature is found in different con¬ 
texts like quantum field theory at finite temperature, on 
non-commutative manifolds and on a curved spacetimes 
(see 46 and references therein). Still, hints of a phase 


transition from a symmetric to a broken phase, in the 
approximation of large size of the group manifold, were 
found. 

Note that progress towards the a better characterisation 
of the phase diagram and of phase transitions in tensor 
models has also been recently achieved [47||48]. In fact, 
using a similar mode integration alongside double scaling 
limit techniques, the nonperturbative analysis of quartic 
tensor models has been performed, with evidence of spon¬ 
taneous symmetry breaking mechanism similar to that 
suggested in [45]. The models considered (quartic tensor 
models with trivial kinetic term) as well as the techniques 
employed (double scaling and intermediate field represen¬ 
tation, allowing to solve quartic tensor models with ma¬ 
trix models methods) are very different from the ones em¬ 
ployed in our present study and in 45 , which are based 



Fig. 1: Colored symmetric interaction terms. 


on generic FRG conventions and concern TGFT models 
with non trivial kinetic kernels. 

In this work, we study a class of TGFT models which 
possess no such additional scale, thus are expected to show 
proper fixed points, so, in a sense, improve on the previous 
analysis. The model we consider is a rank-3 TGFT with 
fields defined on the non-compact manifold R 3 , and en¬ 
dowed with a Laplacian kinetic term. Being a TGFT, this 
model is of interest as a toy model for quantum gravity, 
due to the combinatorics of its Feynman diagrams, but 
even more because it can be seen as a (much) simplified 
version of Lorentzian (T)GFTs for 4d quantum gravity, 
also based on a non-compact group manifold. Thus it can 
be seen as a useful exercise on the way to a renormalisa¬ 
tion analysis of more realistic models, hopefully providing 
useful hints of what to expect for them. 

One certainly generic feature is that the non-compact 
manifold introduces IR divergences, that we properly ad¬ 
dress through a careful definition of a thermodynamic 
limit for TGFTs, in this FRG context. This is an im¬ 
portant technical lesson for later developments. In this 
limit, we recover an autonomous system of /3-functions of 
the coupling constants, and we can then identify the UV 
and IR fixed points of the RG flow. We also find evidence 
for a phase transition (in the continuum limit) from a sym¬ 
metric to a broken (or condensed) phase (which would be 
consistent with the “geometrogenesis” scenario, if we had 
a full geometric interpretation for the simple TGFT model 
we are considering). 


The model. — We consider a rank-3 TGFT defined 
endowed with a specific </> 4 interaction called “mel- 
shown in Fig0 In general, rank-3 melonic 
interactions correspond to peculiar triangulations of the 
3-sphere and are the most dominant objects in the large 


over 

onic ! 


34 


cut-off N limit 23 , 36} 37| (in both simple tensor models 
and topological GFTs, but this result is expected to ex¬ 
tend to a wider class of models). Written in momentum 
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spac^J the classical action of the model reads: 

StyAl = d P 0123(51^ +/ i ) < ^ 123 (!) 

+ ^ J g dpdp' </>12301'23^1'2'3'012'3' + Sy 111 ! 1 , 2 , 3 }| , 


where we used the notation (/> 12 3 = 4 >{PhP 2 iP 3 ) for the 
field modes and “sym” indicates that we include all the 
interactions obtained by symmetrisation over the color la¬ 
bels (see FigQ- The kinetic term is defined by a sum 
of Laplacians acting on the field indices and a mass term 
with coupling p. It is immediate to see that the action is 
built using generalised traces over field indices convoluted 
with the kinetic and interaction kernels and, once expo¬ 
nentiated, defines a quantum theory through a Gaussian 
field measure of covariance (X) s Ps + /r) _1 . 


FRG equations for tensorial models. — The Func¬ 
tional Renormalisation Group approach 38-41 rephrases 
the problem of integrating out the high modes of a theory, 
as one of solving a differential equation, the FRG equa¬ 
tion. Being non-perturbative in nature (with respect to 
any expansion in the interaction coupling constants), the 
FRG allows in principle to deal with the full set of quan¬ 
tum fluctuations of the model and to study its critical 
behavior. 

The implementation of the FRG method in TGFT [45] 
follows closely the usual one [38p9pTl , with special atten¬ 
tion payed to the fact that we are dealing with a convo¬ 
lution of tensors and thus with peculiar non-local interac¬ 
tions. We start by decoupling the field modes as typical in 
the Wilsonian approach to renormalisation, by adding to 
the action a mass-like regulator term A S k = Tr(</> ■ R k ■ <fi) 
depending on the IR cut-off fc, that splits the modes into 
high modes (|p| > k) and low modes (|p| < k). The scale 
dependent quantum theory will then be defined through a 
partition function where high modes are integrated out: 


Z k [J,J\ 


J dcjidc/) e _s ^fo ASfc [^l +ri> ( J ^) +Tr (‘ 7 ' ?i ), 


( 2 ) 


where J is a complex tensor playing the role of a source 
and Tr(J(/>) := J R3 Ji 23 ^i 23 - After a Legendre transform, 
we identify a scale dependent effective action which en¬ 
codes the full information about the quantum theory: 

T k [p,p] = sup{Tr(J^) +Tr(Jv?) - W k [J, J} - AS*,^,^]}, 
j,j 

( 3 ) 

where ip = {<j>) and W k [J, J] = log Z k [J, J]. The term 
A S k is also chosen to be compatible with the choice of 
initial conditions for the FRG differential equation, en¬ 
coding the scaling of effective action to the bare one in 

1 We adopt the standard QFT terminology for field modes, even 

though no spacetime interpretation is associated to the domain of 

the fields, and thus no standard physical interpretation should be 

associated to their modes. The same remark applies to our use of 

the terms ‘UV’ and ‘IR’ throughout the article. 


the UV: T k [<p,<p] —> SVp,p], where A plays the role of 

k — yA. 

a UV cut-off. Introducing the logarithmic scale t = log k 
( 2 ) 

and T k 1 := 5T k /6pSp, the Wetterich equation for tenso¬ 
rial GFT models has the form [45] : 

d t T k = Tr (dtR k ■ [rf + R ^- 1 ). (4) 


Equation Q is fully non-perturbative and exact, and en¬ 
codes (formally) all the information about quantum fluc¬ 
tuation in a typical one-loop form. 

Before moving on to the solution of this equation and 
to the study of the critical points, let us anticipate an 
important technical fact, which we will have to deal with 
in the following. Despite the presence of momentum cut¬ 
offs, evaluating ([4]) requires, as in the ordinary scalar field 
theory, an infinite volume regularisation. In the local field 
theory case, infinite volume divergences are cured by pass¬ 
ing to constants field modes or taking a thermodynamic 
limit 38,39 . In the case of compact groups as worked 
the issue does not arise, as the group vol- 


45 


out in 

ume is finite. The presence of a finite radius results in 
a non-autonomous system of equations where the scale k 
appears explicitly. The existence of phase transition can 
be inferred only in the limit of infinite radius 


46 . The 


present situation differs from both cases, namely the local 
field theory and compact TGFTs in the limit of infinite 
volume. Because of the crucial non-local properties of the 
interactions, the use of constant field modes is misleading 
(indeed (j) 4 terms in the TGFT case do not have the same 
combinatorics, and this cannot be neglected) and, in the 
compact group case, one must understand how to best 
perform, both from a conceptual and practical point of 
view, an infinite radius limit in the equations. In our case, 
we then find wiser to perform a thermodynamic limit. 


Truncation scheme. — In order to be able to per¬ 
form practical computations, we need to adopt a trunca¬ 
tion scheme for the effective action. 

Of course performing a truncation means loosing the 
exact nature of the Wetterich equation. Generally, this 
also generates a singularity of the flow that splits the space 
of couplings in disconnected regions. In a neighbourhood 
of the singularity, we cannot trust the computations. Since 
one usually is interested in the free theory around which 
the perturbative expansion makes sense, we will discuss 
only the region connected to the origin in the space of 
couplings. 

We choose to truncate Tfc to the quadratic term in the 
derivative of the fields and to order four in the fields, thus 
obtaining a form similar to the action itself: 


r fc W, ¥>] = J ^ dp ¥>123 Y,P 2 s+ ^ 123 ( 5 ) 

f dpdp' ^123^1'23^1'2'3'^12'3' + s y m ( 1 > 2 > 3 i • 

Z J K<3 L L J . 


We can already see that, in this way, the UV initial con¬ 
dition on the flow is satisfied. 
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From ([5]), the 2-point 1PI Green function expresses as 
ll 2) (p,p') = (Zk'EsPs + Pk)S{p- p’) + Fk(p,p'), where 


Fk{ p,p') = A fe 


dqi ‘Pqip' 2 p' 3 ‘P qiP 2 P J(Pi -Pi) 


+ j dq 2 dq 3 ip p > iq 2 q 3 ip piq 2 q 3 6 (p 2 - p' 2 ) 6 {p 3 - p' 3 ) 


sym 


{ 1 , 2 , 3 } 


( 6 ) 


The regulator function is chosen as [49]: R k { p, p’) = 5(p— 
p’)Z fc (fc 2 - Yf s =iP^)d{k 2 - Ss =1 Ps), where 9 stands for 
the Heaviside step function. This is a standard choice and 
it satisfies all basic requirements, namely: Rk=o = 0, Vp, 
so that Zk=o[J, J] = Z[J, J ]; Rk =a oc A 2 , Vp s.t. |p| < k, 
to approximately freeze the propagation of modes with 
norm smaller than k; i?fc(|p| > k) = 0, so that high modes 
are unaffected by the regulator. In addition, this choice is 
particularly interesting in our framework because its func¬ 
tional properties allow the analytic evaluation of spectral 
sums. 

If we act on the regulator with the derivative with re¬ 
spect to the logarithmic scale, we find dtR k = d(k 2 — 
E s p 2 )[d t Z k {k 2 - E s p 2 ) + 2k 2 Z k \5(p - p') and the 6 {k 2 - 
E s p 2 )-term so generated simply cancels out. 

Expanding the Wetterich equation, it seems natural to 
choose an expansion in powers of (<p<p), which we perform 
up to the third order, and discarding the vacuum terms, 
to obtain our final truncated functional equation, from 
which we read out the differential equations for the beta 
functions of the theory. 


Thermodynamic limit. — In order to regularise vol¬ 
ume divergences, we perform a lattice regularisation in the 
p-space, which follows from a compactification in the di¬ 
rect space, according to the conventions of [50]. Defining 
the model Q over a lattice D* = [^Z] 3 = [-Z] 3 := [ZZ] 3 , 
of spacing T 3 proportional to the volume of the direct 
space, the Fourier transform becomes a Fourier series 
and, for any function /(p), we have f D , dpi /(p) = 
^ S{ Pi }eu* /(p). We define the delta distribution in D* 
as: 60 * (p, q) = <5 P , q /Z 3 , with d P:q , the Kronecker delta. 
As a result, we have: 8 d* (p,p) = ^p.pA 3 = 1A 3 - 

Using this regularisation prescription, the effective ac¬ 
tion of the model reads: 

T k [ip^-,l] = l 3 ^2 V 123 (Zk^pl + Pk)pi23 (7) 

pS-D* s 

V>i23‘Pi'23 t Pi'2'3'Pi2'3' +sym{l,2,3| , 

where <p(p) = f D [dXi]f =1 e~‘^i piXi ip(x 1 ,x 2 ,x 3 ). In the 
end, the continuous description will be recovered in the 
thermodynamic limit l —> 0. 

The dependence of the system on the volume of the di¬ 
rect space is now explicit, and we can tune this dependence 
in order to consistently remove all the divergences, and be 
left with he physical /3-functions. 



/3-functions. — The IR regularisation of the system of 
/3-functions is direct: we need to extract from the coupling 
constants an explicit dependence on the volume of the di¬ 
rect space, in addition to their scaring with the momentum 
cut-off. After a lengthy but straightforward calculation, 
the set of /3-functions computed with the prescription in¬ 
troduced in the previous section read^J 


P(Z k ) = 


Ai. 


(■ Z k k 2 


Pk ) 2 _ 

. k 2 k 
+dtZk ( tt— + 4 — 


2 Z k 2— + 7T -pr 


r- 


P{Pk) = -3 




{ZkW + pk ) 2 2Zh 
/ 7T k 4 4 fc 3 

+dtZk \ 2 ¥ + 3T 


.k 3 


P 


/3(A k ) = 


2A l 


( Z k k 2 


Pk ) 3 . 
7 r k 4 


2 Zk 


k 4 

¥ 


10 - 


2 k‘ 


„ ^ .. 20 fc 3 

+d t Z k [ 2p+y T 


2 k‘ 


( 8 ) 


To make sense of it in the infinite volume limit, we use the 
ansatz: 


= Z k Rk~x, pk = PkZkRk 2 -*, A k = \kZ\l¥° , 

(9) 

where [Z k \ = [p k \ = [X k ] = 0, [ip\ = -f and £ + er = 4. 
These dimensions are fixed by requiring that [F k ] = 0 with 
[p] = 1. Now, extracting the dimensionless /3-functions, 
one gets: 


_ A k lt k ° 

1 Ik — Pxfc4-2.x( 1+7 j fc )2 


n k k 2 

2[2 j+ 7t ¥ 


P(Pk) = - 


3X^1^^ 


/3(A k ) = 


J3xfe6-3x(l+ Pfc )2 

1 h.3\ 1 

21 n ¥ + 2 T 

2A l^k” 

l 3 xk e ~ 3 x(l+jL k ) 3 


/ \i 4 k k 2 

\Vk ~ X) ( 4 y +7r -[2 


X, 

( 7r k 4 4 k 3 

(%-x)^f + 3T 

- PkPk - (2 - x)Pk , 


. 7T k 4 20 fc 3 2 

(%-x)( 2 ¥ + yy +2fc 


+ 2 ( TT^- +10y +2 ¥ 


- 2p k x k - crAfc . (10) 


The system (101 of /3-functions is non-autonomous in the 


IR cut-off k as long as the parameter l is kept finite. This 
feature is due to the peculiar combinatorics of the ten- 
sorial vertices which span the 1PI 2-point functions with 
different volume contributions. One way to realise this 
is by noting the unusual delta distributions in F in ([6]). 
From (10), we see two different systems arising in the UV 


and IR cut-off limits, coming from different leading terms. 
The fact that the set of /3-functions of a TGFT over a 
compact group manifold is non-autonomous is consistent 


2 We drop here the symbol 1 irri/_ ^q to simplify the notation. 
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with the analysis of standard field theories on compact 


(and curved) manifolds 46 


In order to make sense of the non-compact limit, we 
solve the system in the variables £ and x by requiring 
that the highest volume contribution is regularised and 
all the sub-leading infinities are sent to zero. We have: 
r £ - 2y - 2 = 0 

< which yields v = 0, £ = 2, a = 2. The 

l£-3x-2 = 0 J S 

resulting system of differential equations for the theory is: 


Vk 

< / 3 ( 74 ) 
/3(A k ) 


+ 2 ) 


(1 + Mfe ) 2 

/Vk . 0 X _ 0 — 

: -(TT^ ( T + 2) -’ w ‘*- 2, ‘‘ 

v2 

■ ( , — x 3 (Vk + 4) — 2?7fcA k — 2Afe 

(! +/R-) 3 


( 11 ) 


which is the starting point of our computation of the RG 
flow. As expected, absent any remaining fixed external 
scale, the system is now autonomous. 

The RG flow. — Proceeding with the standard anal¬ 
ysis, we first determine the fixed points and then study 
the linearised system around them to determine the crit¬ 
ical exponents of the model. From the non-linear nature 
of the /3-functions, we have a singularity at JL k = — 1 and 
Afc = (1 + JT k ) 2 /n. In a neighbourhood of those singu¬ 
larities, we do not trust the linear approximation and, 
being interested mainly in the sector of the theory con¬ 
nected with the Gaussian fixed point (i.e. to the pertur¬ 
bative regime of the theory), we will not study the flow 
around points beyond the singularities. By numerical eval¬ 
uation, we find a Gaussian fixed point (GFP) and three 
non-Gaussian (NGFP) fixed points in the plane (ji k , A*,). 
We discard one of them because it lies beyond the singu¬ 
larity. The others correspond to 



Pi = (8.619, -47.049), P 2 = 10 _1 (-6.518,0.096) . 

The stability matrix at the GFP has an eigenvalue with al¬ 
gebraic multiplicity 2 corresponding to the canonical scal¬ 
ing dimensions of the couplings: 8^ 2 = —2, but one sin¬ 
gle eigenvector v = (1,0), thus, considering that all the 
trajectories flow into the origin, the GFP must have a 
marginal direction in the UV. In a neighborhood of the 
non-Gaussian fixed points, we have: 



- 0.351 for v)j_ ~ 

10" 1 (0.65, -9.98), 

(12) 

0 l _ r 

- -2.548 for v]_ 

~ 10 _1 (-6.88,7.26), 

(13) 

e\r- 

- 10.066 for ' 

- 10 -1 (9.996,-0.269), 

(14) 

e 2 _r- 

- -1.988 for v^_ 

~ 10 -1 (9.987,0.506). 

(15) 


The flow of the couplings between the two NGFPs and 
the Gaussian one are plotted in Figj2] The origin is a UV 
sink for the flow; hence, the model is asymptotically free. 
As mentioned before, the absence of a second eigenvector 


Fig. 2: RG flow of the model - The red and blue lines represent, 
respectively, the zeros of /3{jl k ) and /3(Afc), the brown arrows 
are the eigenperturbations of the non-Gaussian fixed points (in 
black), while the green ones those of the Gaussian fixed point 
(in red). Arrows point in the UV direction. The thick black 
line is the singularity of the flow. 

for the stability matrix around the GFP requires an ap¬ 
proximation beyond the linear order and is a signal of the 
presence of a marginal perturbation. By close inspection 
of the plots, confirmed by direct integration at second or¬ 
der of the system of /3-functions, which can be performed 
for generic numerical constants/initial conditions, we infer 
that the behaviour of this direction is still UV attractive, 
i.e. that it corresponds to a marginally relevant direction. 

Both the non-Gaussian fixed points have one relevant 
and one irrelevant direction. They are also characterised 
by the so-called “large river effect”. This effect shows a 
splitting of the space of coupling in two regions not con- 
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nected by any RG trajectory. Thus, the irrelevant di¬ 
rection for the NGFP match the properties of a critical 
surface and suggests the presence of phase transitions in 
the model. In the A > 0 plane, the flow is similar to 
the one of standard local scalar Held theory on R 3 in a 
neighbourhood of the Wilson-Fisher fixed point. That is: 
above the critical surface, the IR limit of the RG trajec¬ 
tories brings the theory in a region where both JL k and Xk 
are positive, while below the irrelevant eigendirection for 
P 2 , the mass parameter is driven to be negative in the IR, 
indicating a spontaneous symmetry breaking mechanism 
(in the different but related context of tensor models, such 
a mechanism has been also found in 48]). In the sector 
A < 0, the situation is rather peculiar. We might infer 
that Pi has the same properties just discussed but re¬ 
versed with respect to the critical surface. The symmetric 
phase where ~p k and A& have the same sign in the IR is 
below the irrelevant direction of the fixed point, while the 
broken phase lies above it. In this sense, we have a phase 
transition also crossing the surface A = 0, but this is not 
an irrelevant direction for any NGFP. This feature sug¬ 
gests that, in this case, we may have a first order phase 
transition. Nevertheless, we must remember that the sec¬ 
tor A < 0 generates theories with a non-stable coupling, 
which is generally not considered in a field theory context. 
This sector must therefore be analyzed under a different 
parametrisation, if we want to shed more light on it. 


In a GFT model with additional geometric data, and a 
proper simplicial gravity interpretation, a broken or con¬ 
densate phase could be interpreted as a continuum geo¬ 
metric phase 19 21 , and would support a geometroge- 


nesis scenario for the emergence of continuum spacetime 
and geometry from these GFT models. The model under 
consideration would therefore need to be enriched with 
such additional data to be more than an indirect support 
for such a scenario. Also in our model, in any case, a 
proper study of the broken phase, involving a change in 
parametrisation for the effective potential and a detailed 
study of the theory around the new ground state, solving 
the classical equation of motion of the model, in a saddle 
point approximation, would be needed to confirm conclu¬ 
sively the existence of a phase transition as envisaged. 


* * * 
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